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ABSTRACT 


A description  of  the  computer  program  used  to  measure  the  sound 
power  output  of  a sound  source  in  a reverberation  room  is  presented. 

The  program  controls  the  operation  of  a real-time  analyzer  and  a micro- 
phone multiplexer.  The  sound  power  is  computed  according  to  ANSI  SI. 21. 
The  program  is  written  in  FORTRAN  V and  requires  approximately  37,000 
eight-bit  bytes  of  core  memory.  Flow  charts,  source  listings,  and  sample 
printouts  are  included. 
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I.  INTRODUCTION 


This  report  describes  the  computer  program  used  to  measure  the  sound 
power  output  of  a sound  source  in  a reverberation  room  according  to  the 
American  National  Standard  Methods  for  the  Determination  of  Sound  Power 
Levels  of  Small  Sources  in  Reverberation  Rooms,  SI. 21-1972.  The  program 
is  written  to  control  the  instrumentation  system  used  for  the  NBS  rever- 
beration room.  It  could,  however,  be  suitable  for  use  in  any  laboratory 
with  a similar  instrumentation  system. 

II.  INSTRUMENTATION 

The  program  assumes  that  two  elements  of  the  instrumentation  system 
are  under  computer  control.  These  elements  are  a real-time  analyzer  and 
a microphone  multiplexer.  The  real-time  analyzer  filters  the  input 
signal  through  30  contiguous  1/3-octave  band-pass  filter  channels, 
samples  each  channel  simultaneously  for  an  integration  period,  determines 
the  rms  level  for  each  channel,  and  transmits  the  results  to  the  computer. 
For  the  NBS  system,  the  time  required  to  transmit  the  levels  in  all  the 
1/3-octave  bands  to  the  computer  is  approximately  1 ms.  The  real-time 
analyzer  integration  time,  under  computer  control,  is  entered  by  the 
operator  from  a CRT  terminal. 

The  computer  controls  the  channel  selection  of  a microphone  multi- 
plexer. In  the  NBS  system,  the  multiplexer  has  16  channels,  of  which  12 
are  used  with  the  12-microphone  array  installed  in  the  reverberation  room. 

III.  PROGRAM  SEQUENCE 

The  execution  of  the  sound  power  measurement  is  summarized  in  the 
abbreviated  flowchart  given  in  Appendix  A.  At  the  beginning  of  the  program, 
the  table  of  reverberation  times  for  the  room  are  entered.  The  operator 
has  the  option  of  entering  the  data  from  the  CRT  terminal  keyboard  or  from 
the  disc  file  output  of  program  REVERB. 


1 


The  microphone  array  calibrations  are  then  acquired.  These  are 
necessary  in  order  to  convert  the  real-time  analyzer  data  from  voltage 
levels,  in  dB  re  1 mV  to  sound  pressure  levels,  in  dB  re  20  yPa.  The 
calibrations  for  each  channel,  in  the  form  of  corrections  to  be  added  to 
the  RTA  data,  may  be  determined  in  advance  and  entered  from  the  CRT 
terminal  keyboard,  or  they  may  be  measured  at  this  point  in  the  program 
by  relating  the  real-time  analyzer  output  to  a reference  signal,  such  as 
an  acoustic  calibrator. 

The  frequency  response  corrections  for  the  microphones  are  then 
entered  from  either  a disc  file  or  from  the  CRT  keyboard.  These  correc- 
tions may  include  the  frequency  response  characteristics  of  the  micro- 
phones and  the  responses  of  the  1/3-octave  frequency  channels  of  the 
real-time  analyzer,  relative  to  the  1000  Hz  band. 

The  atmospheric  pressure,  temperature  and  humidity  (which  the 
program  computes  from  the  wet  and  dry  bulb  temperatures),  the  vane  speed 
and  the  integration  time  are  now  entered.  The  ambient  noise  is  then 
measured,  if  desired,  followed  by  a sound  pressure  level  measurement  at 
each  microphone  with  the  sound  source  operating. 

The  sound  power  levels  are  computed  and  printed  out.  The  operator 
has  the  option  of  either  proceeding  immediately  with  another  measurement 
using  the  same  reverberation  time  data  and  the  same  system  gain  settings, 
or  terminating  the  program. 

IV.  COMPUTATIONS 

The  following  quantities  are  computed: 

1 . Sound  pressure  level  in  the  kth  1/3-octave  band  at  microphone  j. 

Ljk  = Xjk  - Ak  + Bj 


2 


where  X ^ is  the  reading  in  the  kth  1/3-octave  band  at  microphone  j, 
is  the  frequency  response  correction  for  the  kth  1/3-octave  band  for  both 
the  microphone*  and  analyzer,  and  fL  is  the  level  correction  at  1000  Hz  for 
microphone  j . 

The  X.^  are  printed  in  the  columns  labeled  MIC  1 to  MIC  12  in  Printout 
1 of  Appendix  B.  The  A^  are  printed  in  the  column  labeled  FREQUENCY 
RESPONSE  CORRECTION.  The  B_.  are  printed  in  the  column  labeled  CORRECTION 
TO  RTA  DATA.  The  L ^ are  printed  in  the  columns  labeled  MIC  1 to  MIC  12 
in  Printout  2 of  Appendix  B. 

The  uncertainty  in  the  sound  pressure  levels  is  dependent  upon  the 
sampling  characteristics  of  the  real-time  analyzer,  which,  in  the  NBS  sys- 
tem, has  a resolution  of  0.25  dB.  The  levels  are  printed  to  the  nearest 
0.01  dB  in  order  to  aid  in  tracing  the  computations  performed  on  the  data, 
rather  than  to  represent  the  uncertainty  of  the  data. 

2.  Standard  deviation  of  Lj^  (according  to  ANSI  SI. 21). 


a . = 


N. 


m 


-|l/2 


(Nrc-l)'1  J]  O-jk" 
j = l 


^k> ) ' 


where  Nm  is  the  number  of  microphones, 

% 

<Lk>  - ir  Z Ljk 

m -j.i 

*The  program  can  be  modified  to  permit  the  input  of  a separate  frequency 
response  for  each  microphone  channel  to  allow  for  differences  among  the 
microphone  cartridges  used  in  the  array.  As  presently  written  it  is 
assumed  that  all  the  microphones  have  the  same  frequency  response 
corrections . 
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The  are  printed  in  the  column  labeled  STD.  DEV.  in  Printout  4 of 
Appendix  B. 

3.  Mean  sound  pressure  level  in  the  kth  1/3-octave  band. 


The  are  printed  in  the  column  labeled  SOUND  PRESSURE  in  Printout 
4 of  Appendix  B. 

4.  Sound  power  level  in  the  kth  1/3-octave  band. 


Pk  = Lk  - 10  log10  (Tk/T0)  + 10  log1()  (V/V0) 

+ 10  log10  (1  + SAk/8V)  - 10  log  j q (B/1000)  - 14 

where  T^  is  reverberation  time  in  the  kth  1/3-octave  band  in  seconds,  TQ 


wavelength  in  meters  at  the  center  frequency  of  the  kth  1/3-octave  band,  S 
is  the  total  surface  area  of  the  room  in  square  meters,  and  B is  the  baro- 
metric pressure  in  mbars. 

The  P^  are  printed  in  the  column  labeled  1/3  OCTAVE  BAND  SOUND  POWER 
in  Printout  4 of  Appendix  B. 


Lk  = 10  log 
K 10 


is  1 s , V is  the  volume  of  the  room  in  cubic  meters,  V is  1 m3 , A.  is  the 
* ’ok 
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V.  SOFTWARE  DETAILS 


The  programming  language  used  is  the  Interdata  FORTRAN  V Level  1 
Software  System.*  It  is  a superset  of  ANSI  Standard  X3. 9-1966  FORTRAN  and 
should  be  compatible  or  nearly  compatible  with  the  FORTRAN  V used  on  other 
computer  systems.  Users  of  this  program  on  other  systems  should  check  the 
compatibility  of  FORMAT  statements  using  the  A Format  and  INTEGER*  2 
statements  which  define  the  storage  size  of  integers  to  be  two  8-bit  bytes. 

The  program  for  measuring  sound  power  consists  of  a main  program, 
SPOWER,  written  in  FORTRAN  V,  several  FORTRAN  subroutines,  and  several 
assembly  language  subroutines  written  for  use  on  an  Interdata  Model  70 
minicomputer.  These  elements  are  listed  in  Table  1.  The  corresponding 
flow  charts,  sample  printouts,  and  listings  are  given  in  Appendices  B 
and  C. 

In  the  NBS  system,  the  complete  task  occupies  37,066  eight-bit  bytes 
of  memory. 

The  flow  charts,  given  in  the  appendices,  conform  to  the  American 
National  Standard  Flowchart  Symbols  and  Their  Usage  in  Information 
Processing,  X3. 5-1970.  In  these  flow  charts,  the  small  numbers  enclosed 
in  parentheses  refer  to  FORTRAN  statement  labels  in  the  program  listings. 


*The  commercial  computer  products  utilized  are  identified  in  order  to 
adequately  describe  the  program  elements  discussed  in  this  report.  In  no 
case  does  such  identification  imply  recommendation  or  endorsement  by  the 
National  Bureau  of  Standards,  nor  does  it  imply  that  these  products  are 
necessarily  the  best  available  for  the  purpose. 
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In  the  program  listings,  the  first  three  statements  of  every  program 
or  subprogram  are  control  statements  used  to  insert  the  program  name,  in 
the  form  of  a binary  label,  into  the  compiled  object  code. 

VI.  INPUT/OUTPUT  LOGICAL  UNITS 

The  following  logical  units  are  used  for  input/output  operations: 

LU1  - File  containing  a table  of  the  reverberation  time.  This  file  is 
generated  by  the  program  REVERB. 

LU2  - Output  file  for  storing  the  computed  sound  power  levels.  This  file 
is  rewound  at  the  start  of  the  program  and  the  results  of  each 
sound  power  measurement  are  successively  stored  on  it. 

LU3  - Hard  copy  printout  device.  Should  be  assigned  to  a 132-character 
line  printer.  During  execution  of  the  program,  the  operator  may 
choose  the  printout  device  logical  unit  to  obtain  either  a hard 
copy  printout  or  a view  of  the  results  on  the  CRT  terminal. 

LU4  - File  containing  frequency  response  corrections  of  the  microphones. 
These  corrections  consist  of  a set  of  30  numbers,  one  per  band, 
representing  the  difference  in  total  system  sensitivity  between 
each  band  and  band  30  (1000  Hz).  The  microphone  frequency  response 
and  the  differences  in  gain  between  the  real-time  analyzer's  multi- 
filter channels  are  to  be  incorporated  in  these  corrections. 

LU5  - Operator  communication  device.  Should  be  assigned  to  a CRT  terminal. 
Through  this  device  the  operator  controls  the  course  of  the  measure- 
ment, making  decisions  that  are  called  for  and  entering  certain 
measurement  parameters.  These  parameters  include: 

1.  The  reverberation  times,  if  not  read  from  a disc  file. 

2.  The  microphone  calibration  data,  if  not  obtained  during  the 
measurement  with  a reference  signal. 
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3.  The  frequency  response  corrections,  if  not  read  from  a disc  file. 

4.  The  temperature,  relative  humidity,  and  atmospheric  pressure. 

5.  The  real-time  analyzer  integration  time. 

Transmission  of  data  to  and  from  the  real-time  analyzer  and  the 
microphone  multiplexer  is  not  handled  through  logical  units,  but  through 
separate  subroutines  written  in  the  computer's  assembly  language.  These 
subroutines  are  included  in  Table  1. 
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TABLE  1.  SOUND  POWER  PROGRAM  ELEMENTS 
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*Details  of  these  routines  are  given  in  T.  Bartel,  "Interactive  Computer  Program  for  the  Determination  of 
Reverberation  time,"  NBSIR  77-1383,  National  Bureau  of  Standards,  Washington,  D.C.  20234,  December  1977. 


TABLE  1.  (Continued) 
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i 
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♦Details  of  these  routines  are  given  in  T.  Bartel,  "Interactive  Computer  Program  for  the  Determination  of 
Reverberation  Time,"  NBSIR  77-1383,  National  Bureau  of  Standards,  Washington,  D.C.  20234,  December  1977. 
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Abbreviated  Flow  Chart  of  Sound  Power  Program 
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ABBREVIATED  FLOW  CHART:  SPOWER 
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ABBREVIATED  SPOWER-2 
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ABBREVIATED  SPOWER-3 
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APPENDIX  B 


Program  SPOWER  Flow  Chart,  Terminal  Messages,  Printouts  and  Listings 


14 


FLOW  CHART:  SPOWER 
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SPOWER-2 
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SPOWER-3 
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SPOWER-4 
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SPOWER-5 
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SPOWER-6 
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SPOWER-7 
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SPOWER-8 
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SPOWER-9 
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SP0WER-10 
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SPOWER-11 
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SPOWER-12 
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SPOWER-13 


CRT  TERMINAL  MESSAGES:  PROGRAM  SPOWER 
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CRT  TERMINAL  MESSAGES:  PROGRAM  SPOWER  (page  2) 


29 


CRT  TERMINAL  MESSAGES:  PROGRAM  SPOWER  (page  3) 
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CRT  TERMINAL  MESSAGES:  PROGRAM  SPOWER  (page  4) 
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CRT  TERMINAL  MESSAGES:  PROGRAM  SPOWER  (page  5) 
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PRINTOUT  2 - PROGRAM  SPOWER 
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1/3  OCTAVE  BAND  SOUND  POUER  (DB) 
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BAND  FREQ.  SOUND 

NO . PGUEP 


PRINTOUT  4 - PROGRAM  SPOWER 


SOUND  PQUER 


1/3  OCTAVE  OCTAVE 


BAND 

BAND 

BAND 

REVERB 

STD. 

SOUND 

SOUND 

SOUND 

NO. 

FREQ. 

TIME 

DEV. 

PRESSURE 

POWER 

POWER 

(HZ) 

(SEC) 

(DB) 

(DB) 

(DB) 

(DB) 

14 

25 

22.11 

3.65 

56.  17 

58.88 

15 

32 

17.82 

4.  14 

58.70 

61.75 

67.35 

lb 

40 

14.85 

5.41 

61.64 

65.01 

17 

50 

9.51 

3.01 

61.77 

66 . 65 

18 

63 

4.  14 

1.32 

62.56 

70.66 

76.07 

19 

80 

4.57 

2.04 

66.50 

73.84 

2U 

100 

3.02 

1.37 

67.42 

76.29 

21 

125 

3.71 

0.71 

68.66 

76.40 

31.59 

22 

160 

4.51 

0.30 

70 . 96 

77 . 64 

23 

200 

5. 19 

0.82 

72.97 

78.88 

24 

250 

6.31 

0.44 

74.86 

79.79 

84.  16 

25 

315 

6.80 

0.36 

74.95 

79.44 

26 

400 

6.97 

0.36 

74.83 

79.  13 

27 

500 

6.76 

0.40 

74.93 

79.28 

84.31 

28 

630 

6. 18 

0.41 

75 . 45 

80. 14 

29 

800 

5.91 

0.37 

78.31 

83.  14 

30 

1000 

5.94 

0.42 

79.75 

84.53 

S9.24 

31 

1250 

5.92 

0.36 

80.67 

85.43 

32 

1600 

5.66 

0.37 

79.97 

84.90 

33 

2000 

5.  19 

0.39 

79.35 

84 . 63 

89. 0S 

34 

2580 

4.63 

0.40 

77.43 

83.20 

35 

3150 

4.01 

0.40 

74.97 

81.35 

36 

4000 

3.34 

0.37 

74.39 

81.56 

86.03 

37 

5000 

2.70 

0.39 

72.76 

S0.34 

38 

6300 

2.  12 

0.42 

70.55 

79.67 

39 

8U00 

1.58 

0.49 

67.86 

78.25 

83.10 

40 

10000 

1 . 16 

0.68 

64.78 

76.51 

41 

12500 

0.84 

0.91 

61.53 

74.69 

42 

16000 

0.61 

1.21 

57.50 

72.03 

43 

20000 

1.59 

53.47 

A -LEIGH  TED 

SOUND  POLER 

- 

94.42  DB 

( 2.77 

MILLIWATTS) 

UNWEIGHTED 

SOUND  POWER 

- 

94.80  DB 

( 3.02 

MILLIWATTS) 

A-UE IGHTED 

SOUND  PRESSURE  = 

88.89  DB 

UNWEIGHTED 

SOUND  PRESSURE  = 

89.21  DB 

VANE 

SPEED 

= 7.5  RPM 

NO . 

OF  MICROPHONES  = 12 

SAMPLE  TIME  AT  EACH  MICROPHONE:  30  SAMPLES  X 1.00  SECONDS  = 30.00  SECONDS 


DRY  BULB  TEMPERATURE 
UET  BULB  TEMPERATURE 
RELATIVE  HUMIDITY 
BAROMETRIC  PRESSURE 


21.8  DEG.  C 

16.7  DEG.  C 

59.7  % 

980.8  MILLIBARS 


(71.2  DEG.  F) 
(62.0  DEG.  F) 

(735.7  MM  HG) 


SOUND  SOURCE:  B & K TYPE  4204  REFERENCE  SOUND  SOURCE 


DATE:  3/18/77 

TIME:  7:50  A.M. 
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PROGRAM  SPOWER 


1 $ASSM 

2 SPOUER  PROG  SPOUER  - PROGRAM  TO  MEASURE  SOUND  POUER 

3 SPORT 

4 C 

5 C SUBPROGRAMS  CALLED:  SAMPLE, SPLCAL, DELAY, SUBMUX, MUX, HUMID, BAND, 

6 C MEANLG, MEANAR, RTA 

7 C 

8 DIMENSION  X(50) , RDATA (30, 12) ,SPL(30, 12) , STD (30) ,SPLM(30) ,T60(30) , 

9 1POUER C30).UE IGHT<30), FLOOR (30), SPLCORC 12), RELCAL (12), ARRAY ( 12) , 

10  2S0URCEC 18) -DATE (5) -TIME (5) , DATE0C5) , TIME0 (5) , FRCOR (30) 

1 1 INTEGERS  X0(  1560) -FREQ (30), B,OVFLOU,FF, TX100, LINE ( 101)  , 

12  1DOT, STAR, SPACE, LCARET,  RCARET 

13  DATA  DOT, STAR, SPACE, LCARET, RCARET, FF/X' 2E00’ , X' 2AO0' , X'  2000' , 

14  IX' 3C00' , X' 3E00' , X'  0C00 ’ / 

15  DATA  FREQ/25, 32, 40, 50, 63, 80,  100,  125,  160,200,250,315,400,500, 

16  1630,800, 1000, 1250, 1600,2000,2500,3150,4000,5000,6300,8000, 

17  210000,12500,16000,20000/ 

18  DATA  UE IGHT/-44. 7, -39 . 4, -34.6,-30 .2,-26 .2, -22. 5, -19. 1,-16. I, 

19  1-13.4,- 10, 9, -8. 6, -6. 6, -4. 8, -3. 2,-1 .9, -.8,0. , .6,  1 . , 1 .2,  1 .3, 

20  21 .2, 1., .5,-. 1,-1. 1,-2. 5, -4. 3, -6. 6, -9. 3/ 

21  URITE  (5,  10) 

22  10  FORMAT  (X/X/'THIS  PROGRAM  RUNS  AND  ANALYZES  ONE  SOUND  POUER 

23  1 MEASUREMENT  RE  S1.21'/'UITH  1 SOURCE  LOCATION  AND  UP  TO  12 

24  2 MICROPHONE  POS IT  IONS' /X/X/ 

25  3' THE  FOLLOUING  LOGICAL  UNIT  ASSIGNMENTS  ARE  REQU IRED : ' /X/ 

2b  45X, ' LU 1 — INPUT  FILE  CONTAINING  REVERBERATION  TIME  DATA'/ 

27  55X, ' LU2  — OUTPUT  FILE  FOR  STORING  SOUND  POUER  DATA'/ 

28  65X, ' LU3  — HIGH  SPEED  PRINTER'/ 

29  75X,  ' LU4  — INPUT  FILE  CONTAINING  FREQUENCY  RESPONSE  CORRECTIONS'/ 

30  85X, ' LU5  — CRT  TERMINAL' /X/X) 

31  20  FORMAT  (A4) 

32  REIJIND  2 

33  30  URITE  (5,40) 

34  40  FORMAT  ('DO  YOU  UANT  TO  INPUT  THE  REVERBERATION  TIMES  FROM 

35  1 THE  DISC?') 

36  READ  (5,20)  QS 

37  IF  (QS.EQ.'NO')  GO  TO  100 

38  REUIND  1 

39  READ  (1,20)  C 

40  DO  90  1=1,30 

41  READ  (1,80)  J, J,C,T60(I),C,C, J 

42  80  FORMAT  ( 12, X, 15, X, 2 (F8 . 4, X) , F9 . 5, X, F6 . 3, X, 13) 

43  90  CONTINUE 

44  GO  TO  140 

45  100  DO  130  1=1,30 

46  B = 1 + 13 

47  URITE  (5,  1 10)  B 

48  110  FORMAT  ('ENTER  REVERBERATION  TIME  FOR  BAND  ' , 12, 4X, ' (XX. XX) ' ) 

49  READ  (5,120)  T60(I) 

50  120  FORMAT  (F10.5) 

51  130  CONTINUE 

52  140  DO  160  1=1,15 

53  B = 1 + 13 

54  J=B+ 15 

55  URITE  (5,150)  B, T60 ( I) , J, T60 ( 1+15) 

56  150  FORMAT  ('BAND  ',12,'  — T =' , F6 . 2, 15X, ' BAND  ’,12,'  — T =',F6.2) 

57  160  CONTINUE 

58  URITE  (5,170) 

59  170  FORMAT  (X/'DO  YOU  UANT  TO  RE-ENTER  THESE  DATA?') 

60  READ  (5,20)  QS 
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SPOWER-2 


6 1 
62 
b3 

64 

65 

66 
6 { 

6 a 

69 

70 

71 
74 

74 

~7  C 
< -J 

7 6 


180 

198 

200 

205 

21b 


79 

80 

81  220 


•DO 

84 

85 

86 
8 1 
83 
39 

90 

91 

92 

93 

94 

95 

96 
9 < 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 
109 
1 in 
1 1 1 
1 12 
1 13 
1 14 
1 15 
1 16 
1 17 
1 18 
1 19 
120 


230 


250 


260 

270 

300 

310 


2 b 


330 

400 

•405 


410 

415 


IF  (QS.NE.'NO')  GO  TO  30 
DO  180  1=1,30 
FLOOR ( I) =0. 

DO  ISO  J = 1 , 12 
PDATA ( I , J) =0 . 

SPL ( I, J) =0. 

CONTINUE 
URITE  (5,190) 

FORMAT  CENTER  NO.  OF  MICROPHONES  IN  MICROPHONE  ARRAY 
1 (2  DIGITS)') 

READ  (5,200)  NMIC 
FORMAT  (12) 

RELCAL ( 1 ) =0 
LIRITE  (5,210)  NMIC 

FORMAT  ('ENSURE  THAT  THE  MULTIPLEXER  DIG ISUITCH  IS  SET  TO  ' , I2/X/ 
1 1 DO  YOU  WANT  TO  CALIBRATE  MICROPHONE  CHANNEL  NO.  1'/ 

2' WITH  A REFERENCE  SIGNAL?') 

READ  (5,20)  US 
IF  (QS.NE.'NO')  GO  TO  230 
WRITE  (5,220; 

FORMAT  ('ENTER  THE  CORRECTION  TO  BE  ADDED  TO  THE  RTA  DATA'/ 

1 ' TO  GET  SOUND  PRESSURE  LEVELS  FOR  CHANNEL  1 (XX.XX)') 

READ  (5.120)  SPLCOR ( 1) 

GO  TO  300 
CALL  SUBMUX.C  1 ) 

CALL  SPLCAL (SPLCOP ( 1 ) , 50, X, X0, 1560) 

WRITE  (5,250)  NMIC 

FORMAT  ( ‘ DO  YOU  WANT  TO  CALIBRATE  CHANNELS  2 THROUGH  ',12, 

1'  IN  THE  SAME  MANNER?') 

READ  (5,20)  QS 
IF  (QS.EQ.'NO')  GO  TO  300 
DO  270  I =2, NMIC 
LIRITE  (5,260)  I 

FORMAT  (X/X/' MICROPHONE  CHANNEL  ',12) 

CALL  SUBMUXC I ) 

CALL  SPLCAL (SPLCOR (I), 50, X,X0, 1560) 

RELCAL ( I ) =SPLCOR ( 1 ) -SPLCOR ( I ) 

CONTINUE 
GO  TO  400 

WRITE  (5,310)  NMIC 

FORMAT  (X/'THE  CALIBRATIONS  FOR  CHANNELS  2 THROUGH  ',12, 

1'  RELATIVE  TO  CHANNEL  l'/'MUST  BE  ENTERED'/ 

2' (CALIBRATION  FOR  X = RESPONSE  OF  X - RESPONSE  OF  1)'/X) 

DO  330  I =2, NMIC 
WRITE  (5,320)  I 

FORMAT  ('ENTER  CALIBRATION  FOR  CHANNEL  ',I2,5X, 
l'(+  OR  - XX.XX)') 

READ  (5,120)  RELCAL ( I ) 

SPLCOR ( I ) =SPLCOR ( 1 ) -RELCAL ( I ) 

CONTINUE 
URITE  (5,405) 

FORMAT  (X/X/T7, 'MICROPHONE  CALIBRATION  DATA'/X/X/ 

1' MICROPHONE' ,T15, 'CORRECTION' ,T30, 'CALIBRATION'/ 

2T3, 'CHANNEL' ,T 19, 'TO', T30, 'RELATIVE  TO'/ 

3T3, 'NUMBER' , T16, 'RTA  DATA' , T3 1 ,' CHANNEL  1 ' /X) 

DO  415  1=1, NMIC 

WRITE  (5,410)  I, SPLCOR ( I) , RELCAL ( I) 

FORMAT  (T5, 12, T16, F7 . 2, T32, F6 . 2) 

CONTINUE 
WRITE  (5,170) 
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121 

READ  (5,20)  QS 

122 

IF  CQS.NE.'NO') 

GO  TO  205 

123 

DO  420  1=1,30 

124 

FRCOR ( I) =0. 

125 

420 

CONTINUE 

1 2b 

URITE  (5,425) 

127 

425 

FORMAT  ('DO  YOU 

UANT  TO  CORRECT  FOR  THE  FREQUENCY  RESPONSE 

128 

1 OF  THE'  /'  RTA  MULTIFILTER  AND/OR  MICROPHONE?') 

129 

READ  (5,20)  QS 

130 

IF  CQS.FQ.'NO') 

GO  TO  490 

131 

428 

URITE  (5,430) 

1 32 

430 

FORMAT  ('DO  YOU 

UANT  TO  INPUT  THE  FREQUENCY  RESPONSE 

133 

134 

135 

136 

137 

138 

139 
148 

141 

142 

143 

144 

145 
14b 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 


440 

445 

450 

455 


460 


465 


470 

475 


490 


164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175  525 

176  530 

177  540 

178 

179 

180 


100 

110 


20 


1 CORRECTIONS  FROM  THE  DISC?') 

READ  (5.-20)  QS 
IHDTR=0 

IF  C QS . EQ . ' NO ■ ) GO  TO  440 

INDTR= 1 

REUIND  4 

READ  (4)  FRCOR 

READ  (4)  DATE0,TIME0 

GO  TO  465 

iJPITE  (5,445) 

FORMAT  ('FREQ.  RESPONSE  CORRECTION (BAND  X)  = 

1 RESPONSE (BAND  X)  - RESPONSE (BAND  30) ' /X) 

URITE  (5,455) 

FORMAT  ('ENTER  BAND  NO.,  CORRECTION  (XX,YY.YY)'/ 

1' (ENTER  00,0.  IF  NO  FURTHER  CORRECTIONS  ARE  DESIRED)') 

READ  (5,460)  B,I,C 
FORMAT  ( 1 2 - A 1 , F 1 0 . 5 ) 

IF  (B.EQ.O)  GO  TO  465 

FRCOR (B-  13) =C 

GO  TO  450 

DO  475  1=1,15 

B = 1 + 13 

J=B+15 

URITE  (5,470)  B, FRCOR (I) , J, FRCOR ( 1 + 15) 

FORMAT  (2 ('BAND  ',12,':  CORRECTION  = ',F7.2,9X)) 

CONTINUE 
URITE  (5,170) 

READ  (5,20)  QS 
IF  CQS.NE.'NO')  GO  TO  428 

CALL  H U M I D ( H U MD  T Y , D TE  MR  F , D TE  MP  C , UTE  MR  F , UTE  MP  C , 

1 PRESS  1 , PRESS2, VSOUMD) 

URITE  (5,740) 

READ  (5,  120)  VANE 
URITE  (5,510) 

FORMAT  ('ENTER  THE  RTA  INTEGRATION  TIME  IN  SECONDS  (XX. XX)') 
READ  (5, 120)  T 

IF  ( (T.EQ. . 1) .OR. CT.EQ. .25) .OR. (T.EQ. .5) .OR. (T. EQ. 1 . ) .OR. 

1 ( T. EQ . 2 . ) . OR . ( T . EQ . 4 . ) . OR . ( T . EQ . 8 . ) . OR . ( T . EQ . 16.) . OR . 

2 (T.EQ. 32. ) ) GO  TO  525 
URITE  (5.520) 

FORMAT  ('YOU  HAVE  ENTERED  AN  IMPROPER  INTEGRATION  TIME') 

GO  TO  500 
TX100=T*100 . 

URITE  (5,540) 

FORMAT  ('ENTER  THE  NO.  OF  CONSECUTIVE  SAMPLES  TO  BE  TAKEN ',4X, 

1 ' (2  DIGITS,  MAX=50) ' ) 

READ  (5,2 0 0 ) N S A MP 

IF  (NS AMR . GT. 50)  GO  TO  530 
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181 

545 

C=T*NSAMP 

182 

URITE  (5,550)  NSAMP , T, C 

183 

550 

FORMAT  (X/' TOTAL  INTEGRATION  TIME: 

' . 12, ' SAMPLES  X ' ,F5 

184 

1 ' SECONDS.  = ' , Fb . 2, ' SECONDS ' /X) 

185 

URITE  (5,560) 

lob 

560 

FORMAT  ('DO  YOU  UANT  TO  CONDUCT  A 1 

NOISE  MEASUREMENT?') 

187 

READ  (5,2CDQS 

188 

IF  (QS.EQ.'NO')  GO  TO  650 

189 

CALL  SUBMUX(l) 

190 

URITE  (5,570) 

191 

570 

FORMAT  ('ENTER  "GO"  TO  START  NOISE 

MEASUREMENT' ) 

192 

READ  (5,20)  QS 

193 

CALL  SAMPLE  (0,  NSAMP,  TX1 00, OVFLOU,Xi 

0,  1560) 

194 

DO  600  1=1,30 

195 

B = I + 1 3 

l y b 

CALL  BAND (B, NSAMP, X,X0, 1560) 

197 

CALL  MEANLG  ( NSAMP , X,  FLOOR  (ICC) 

1 98 

600 

CONTINUE 

199 

650 

URITE  '.5,660) 

200 

6 B 0 

FORMAT  ('ENTER  "GO"  TO  START  SOUND 

POLJER  MEASUREMENT' ) 

20  I 

670 

REA  D 1 . 5 , 2 0 ) U S 

202 

DO  710  J= 1 , NMIC 

203 

URITE  (5,680)  J 

204 

680 

FORMAT  ('MICROPHONE  ’,12) 

205 

CALL  SUBMUX(J) 

286 

CALL  DELAYC500) 

207  CALL  SAMPLE (0, NSAMP, TX 100, OVFLOU,  X0, 1560) 

208  IF  (OVFLOU. EQ.0)  GO  TO  700 

209  UR  I TE  ( j , 6 9 0 J 

210  690  FORMAT  ( ''  THE  SYSTEM  GAIN  MUST  BE  RE-ADJUSTED ' ) 

2 1 1 GO  TO  205 

2 1 2 700  DO  710  1 = 1,30 

213  B= 1+13 

214  CALL  BAND (B, HSAMP, X, X0, 1560) 

215  CALL  MEANLG (NSAMP, X, RDATA ( I , J) , C) 

216  710  CONTINUE 

217  UR  I TE  (5,720) 

218  720  FORMAT  (X/X/’ SOUND  POLJER  MEASUREMENT  IS  FINISHED’) 

219  GO  TO  750 

220  730  CALL  HUM  1 D ( HUMDTY,  DTEMPF,  DTE!  1PC , UTEMPF , IJTEMPC,  PRESS  1 , PRESS2 

221  1VSOIJND) 

222  UP  I TE  (5  , 740) 

223  740  FORMAT  CENTER  VANE  SFEED  IN  REV/M  IN  (XX.  X) ' ) 

224  READ  (5,120)  VANE 

225  75U  UR  I TE  (5,800) 

226  800  FORMAT  ( ' IDENTIFY  SOUND  SOURCE  (LIMIT:  1 LINE)') 

227  READ  (5,810)  (SOURCE  Cl), 1 = 1, 1S9 

2 2 8 8 1 0 F 0 R M A T ( 1 8 A 4 ) 

229  UR  1 TE  (5,820) 

230  820  FORMAT  (X/’ ENTER  DATE’) 

231  READ  (5,-830)  (DATE ( I) , I =1,5) 

232  330  FORMAT  (5A4) 

233  URITE  (5,840) 


234  840  FORMAT  (’ENTER  TIME  OF  DAY') 

235  READ  (5,830)  (TIME ( I ) , I = 1 , 5) 

236  VOLUME =42 4. 8 

237  AREA =343. 73 

238  DO  950  1=1,30 

239  DO  900  J = 1 , NM I C 

240  SPL ( I , J ) =RD ATA ( I , J ) +SPLCOR ( J ) -FRCOR ( I ) 
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242 

243 

244 

245 

24b 

247 

24b 

249 

250 

251 

■“ico 

£_U<_ 

253 

254 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

268 

269 

270 

27  1 

272 

2 i 3 

2 i 4 

275 

276 

277 

278 

279 

280 

28  1 

282 

283 

2b4 

285 

286 

287 

288 

289 

29C1 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 
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ARRAY(J) = SPL< I, J) 

900  CONTINUE 

CALL  MEANAR(NMIC. ARRAY, C, STD (I)) 

CALL  MEANLG  ( NMI C,  ARRAY,  SF'LMC  I)  ,C) 

IF  CT60 ( I ) . EQ . 0 . ) GO  TO  950 
U A VE  L G = VS  0 U N D /F  R E Q C I ) 

POUERU)  =SPLM(  I)  + 10*ALOG  10 ( VOLUME *(  1+CAREA*UAVELG)/(8*V0LUME)  ) 
1*1000. /(PRESS2*T60( 1)55-14 
950  CONTINUE 
1000  URITE  (5, 1002) 

1002  FORMAT  ('ENTER  LOGICAL  UNIT  FOR  DESIRED  PRINTOUT  DEVICE') 

READ  (5,  1005)  LU 
1005  FORMAT  (ID 

URITE  (LU,  1010)  (1,1  = 1,12) 

1010  FORMAT  (X/T62, ' RTA  DATA  READOUT' , T121, 'FREQUENCY'/ 

1 ' BAND' ,T16, 'NOISE' , T12 1 , 'RESPONSE'/ 

2T2,  'NO. ' ,T8, 'FREQ. ' ,T16, 'LEVEL' ,T24,9(' MIC  ',11, 3X5 , 

33 ('MIC  ' , 12, 2X) ,T1 21, 'CORRECT ION'  /X) 

DO  1030  1=1,30 
8=1+13 

URITE  (LU,  1020)  8, FREQ ( I), FLOOR ( I) , (RDATA ( I , J) , J= D 12) , 
1FRCORCI) 

1020  FORMAT  (T2,  12, T7, 15, T15, F6 . 2, T23, 12 (F6 . 2, 2X) , T122, F6 . 2) 

1030  CONTINUE 

IF  (INDTR.EQ.0)  URITE  (LU, 1032) 

IF  (LU.EQ.5)  GO  TO  1034 

IF  (IHDTR.EO.l)  URITE  (LU, 1033)  DATE0, TIME0 

1032  FORMAT  (X/X/X) 

1033  FORMAT  (X/T90, ' FREQUENCY  RESPONSE  MEASURED:  ' , 5A4/ 

IT  120, 5A4) 

1034  URITE  (LU, 405) 

DO  1035  I = 1 , NMIC 

URITE  (LU, 410)  I , SPLCOR ( I) , RELCAL ( I ) 

1035  CONTINUE 

URITE  (LU, 1040)  FF 
1040  FORMAT  (A1) 

URITE  (LU, 1 100) 

1100  FORMAT  (X/X/X/T49, 'SOUND  PRESSURE  LEVELS') 

URITE  (LU,  1 105)  (1,1  = 1,12) 

1105  FORMAT  (X/X/X/X/' BAND' , T16, 'NOISE'/ 

1T2, 'NO. ' ,T8, ' FREQ . ' , T16, ' LEVEL' , T24, 9 ( ' MIC  ',11, 3X) . 

23 ('MIC  ' , 12, 2X) /X) 

DO  1110  1=1,30 
8 = I + i 3 

C=FLOOR ( I ) +SPLCOR ( 1 ) -FRCOR ( I ) 

URITE  (LU,  1107)  B, FREQ ( I ) , C, (SPL ( I , J) . J = 1 , 12) 

1107  FORMAT  (T2, 12, T7, I5,T15,F6.2,T23, 12CF6.2, 2X) ) 

1110  CONTINUE 

URITE  (LU, 1040)  FF 
C = 0 . 

DO  1150  1=1,30 

IF  (T60( I) . EQ . 0 . ) GO  TO  1 150 

IF  (POUER ( I ) . GT. C)  C=POUER ( I ) 

1150  CONTINUE 

NSC ALE = 10* ( INTCC/10. )+l) 

IF  (NSCALE.L T.50)  NSCALE=50 
J=NSCALE-50 

URITE  (LU, 1 160)  (I,  I = J, NSCALE . 10) 

1160  FORMAT  (X/X/X/T46, ' 1/3  OCTAVE  BAND  SOUND  POUER  (DB)'/ 

1 X/X/X/T 14,5(13, 17X) , I3/T15, 10(' . ' , 9X) , ' . ' /X) 
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301 

302 

303 

304 

305 
30b 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 
3 J 9 

320 

321 


3 4 3 

324 

325 
328 
327 
348 
3 2 9 


330 

331 

332 


>34 


33b 


338 
3 j>9 

340 

341 

342 

343 

344 

345 
b46 

347 

348 

349 

350 

351 

352 
553 
354 

T cr  c 
jJJ 

356 


359 

360 


DO  1220  1=1,30 
B=I+13 

DO  1 170  J-l,  101 
LINE(J)  =SF'ACE 
1170  CONTINUE 

IF  (T60 ( I ) . NE . 0 . ) GO  TO  1 175 
LJRITE  (LU,  1 172)  B.FREQC  I) 

1172  FORMAT  CI2.T7, 15) 

GO  TO  1220 

1 175  J = IHT (2 . *PQLJER  C I ) +. 5)  -2* (NSC ALE-50)  +1 
IF  CJ.GE.l)  GO  TO  1180 
LI  NEC  T)  =LCARET 
GO  TO  1200 

1 180  IF  CJ.LE.  101)  GO  TO  1 190 
LI NEC  101) =RCAPFT 
GO  TO  1200 
1190  L INE  C J) =STAR 

1200  LJRITE  CLU  . 1210)  B,  FREQ  Cl),  (L  I NE  C J) , J= 1 , 101).  POLER  ( I ) 

1210  FORMAT  C I2.T7, I5.T15, 101A1.T120.F6.2) 

1220  CONTINUE 

LJRITE  CLU,  1230)  FF 

1230  FORMAT  <X/'  BAND’- , T7 , ' FREQ T12 1 . ' SOUND'  / 

1 T2 , ' NO . ' , T 1 2 1 , ' POUER ' /A  1 ) 

URITE  CLU, 1400) 

1400  FORMAT  (T23, 'SOUND  POIJER'/X/X/ 

1T41 , ‘ 1/3  OCTAVE', T54, ' OCTAVE'/ 

2T44, 'BAND' ,T54, 'BAND'/ 

3' BAND ’ , Tib, 'REVERB' ,T25, 'STD. ' ,T34, 'SOUND' , T44. ' SOUND' , T54, 
4'  SOUND'/ 

5T2, 'NO.' ,T7, ' FREQ . ' , T16, ' TIME' . T25, ' DEV. ' , T33, ' PRESSURE' » 
6T44, ' POLJER' , T54, ' POUER'  / 

717 , ' CHZ) ' , Tib, ' CSEC) ' , T25, ' CDB)',T34, ' (DB)’.T44. ' (DB)', 
8T54, ' CDB) '/X) 

DO  1500  1=1,30 
B=I+13 

IF  (T60CI) . NE . 0 . ) GO  TO  1420 

URITE  CLU, 1410)  B , FREQ  Cl), STD  Cl), SPLMC I ) 

1410  FORMAT  CT2, 12, T6, 15. T24.F5 . 2, T33, F6 .2) 

GO  TO  1500 

1420  IF  (INTCB/3.+.9) .EQ. INTCB/3.+.001))  GO  TO  1450 

1430  LJRITE  CLU,  1440)  B,  FREQ  C I ) , T60  C I ) , STD  ( I ) , SPLMC  I ) , POUER  ( I) 

1440  FORMAT  CT2,  12. Tb. 15, T 15, F6 . 2, T24, F5 . 2- T33, Fb . 2. T43, F6 . 2) 

GO  TO  1500 

1450  IF  C cTbOC I- 1 ) . EQ . 0 . ) . OR . CT60C 1 + 1) .EQ.0. ) ) GO  TO  1430 

C = 10  . *ALOG  10  C 10  . **( POUER  C I-D/10. ) + 10 .**( POUER ( I)/10. ) + 

1 10.**CPOUERC  I + D/10. ) ) 

URITE  CLU,  1 460) B,  FREQ C I) ,T60( I) ,STDC I) .SPLMC I) , POUER C I) ,C 
1460  FORMAT  CT2,  12, T6, 15 , T15. F6 . 2, T24, F5 . 2, T33, Fb . 2, T43, Fb . 2, 
1T53.F6.2) 

1500  CONTINUE 
X( 1 ) =0. 

XC2) =0. 

XC3) =0. 

XC4) =0. 

DO  1530  B = 14, 43 
I =B-  1 3 

IF  CT60C I) . EQ . 0 . ) GO  TO  1520 

XC 1) =XC 1)+10.**((POUER<I5+UEIGHT(I))/10.) 

XC2)  =XC2)+10,H°K(POUERCI)/10.) 

1520  XC3)-XC3)+10.**( (SPLMC I)+UEIGHT( I) )/10. ) 
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3 b 1 

X ( 4) =X ( 4) +10 . **(SPLM( I ) / 1 0 . ) 

362  1530 

CONTINUE 

363 

X(5) =10.*ALOG10(X( 1) ) 

364 

X(6) =X( 1)*1 .E-9 

365 

X(7) = 10 . +ALOG lO (X(2) ) 

366 

X ( s ) =X(2) *1 . E ~3 

367 

X ( 9 ) = 10.  *ALOG 10 (X(3 ) ) 

bbo 

X ( 1 0 ) = 1 0 . *ALOG 1 0 ( X ( 4) ) 

3fc.y 

WR  I TE  C LU , 1 540 ) ( X ( I ) , I =5, 10 ) 

370  1540 

FORMAT  ( X/X/ ' A - WE  I G H TE  D S 0 U N D 

371 

14X,  ' ( ' , F8 .2, ' MILLIWATTS) ' / 

= ',F6.2, 


DB' 


DB' ,4X, ' ('  ,F8.2, 

DB'  / 

DB'  /X) 

DTEMPC,  DTE  IIP  F 


SAMPLES  X '-F5.2, 


F 4 . 1 ' DEG,.  F)  ' ) 


,F4.  1.'  DEG.  C',7X/(’,F4.  1, 


1=1,18), 


372  2' UNWEIGHTED  SOUND  POWER  = ' ,F6.2, 

373  3'  MILLIWATTS)'/ 

374  4' A- WEIGHTED  SOUND  PRESSURE  = ',F6.2, 

375  5 ‘ UN WE IGHTED  SOUND  PRESSURE  = ',F6.2, 

376  C=T*NSAMP 

377  WRITE  (LU,  1550)  VANE,  NMIC,  NSAMF1 , T,  C, 

378  1550  FORMAT  ('VANE  SPEED  = ',F4.1,'  PPM'/ 

379  1 ' NO . OF  MICROPHONES  = ',12/ 

380  2' SAMPLE  TIME  AT  EACH  MICROPHONE:  ',12, 

381  3'  SECONDS  = ',F6.2,'  SECONDS' /X/ 

382  4' DRY  BULB  TEMPERATURE  = ',F4. 1, ' DEG.  C',7X, '( 

383  IF  (UTEMPC . EQ . - 1000 . ) GO  TO  1570 

384  WRITE  ( L U , 1 5 6 0 ) WTE MP C , WTE MP F 

385  1560  FORMAT  ('WET  BULB  TEMPERATURE  = 

386  1 ' DEG . F ) ' ) 

387  1570  WRITE  (LU,  1580)  HUMDTY, PRESS2, PRESS  1 , (SOURCE ( I 

388  1 (DATE ( I) , 1=1,5), (TIME(I), I=1,3),FF,FF 

389  1580  FORMAT  ('RELATIVE  HUMIDITY  = ',F4.1,'  %'/ 

390  1' BAROMETRIC  PRESSURE  = ',F6.1,'  MILLIBARS  C, 

391  2F5.1,'  MM  HG)'/X/ 

392  2' SOUND  SOURCE:  ' , 18A4/X/' DATE : ',584/' TIME:  ',584/81/81) 

393  WRITE  (5,1700) 

394  170O  FORMAT  ('DO  YOU  WANT  ANOTHER  PRINT  OF  THESE  RESULTS?') 

395  READ  (5,20)  QS 

396  IF  (QS.NE.'NO')  GO  TO  1000 

397  WRITE  (5,1710) 

398  1710  FORMAT  ('DO  YOU  WANT  TO  CORRECT  ANY  TYPING  MISTAKES  IN  THE 

399  1 PRINTOUT?') 

400  READ  (5,20)  QS 

401  IF  (QS.NE.'NO')  GO  TO  730 

402  WRITE  (2)  POWER 

403  WRITE  (5,1800) 

404  1800  FORMAT  ('DO  YOU  WANT  TO  CONDUCT  ANOTHER  SOUND  POWER  MEASUREMENT' 

405  l'USING  THE  SAME  REVERBERATION  TIME  DATA  AND  THE  SAME'/ 

406  2' SYSTEM  GAIN  SETTINGS?') 

407  READ  (5,20)  QS 

408  IF  (QS.EQ.'NO')  GO  TO  1820 

409  WRITE  (5,1810) 

410  1810  FORMAT  ('DO  YOU  WANT  TO  USE  THE  SAME  TEMPERATURE, HUMID  I TY, ' / 

411  1 ' VANE  SPFED,  AND  INTEGRATION  TIME?') 

412  READ  (5,20)  QS 

413  IF  (QS.EQ.'NO')  GO  TO  490 

414  GO  TO  545 

415  1820  WRITE  (5,1830) 

416  1830  FORMAT  (X/X/'END  OF  JOB') 

417  END 
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Flow  Charts,  Terminal  Messages  and  Listings 
for  Subroutines  SPLCAL,  MEANLG,  and  ME ANAR 


FLOW  CHART:  SUBROUTINE  SPLCAL 
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CRT  TERMINAL  MESSAGES:  SUBROUTINE  SPLCAL 
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CRT  TERMINAL  MESSAGES:  SUBROUTINE  SPLCAL  (page  2) 
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SUBROUTINE  SPLCAL 


1 $assm 

SPLCAL  PROG  SPLCAL  - SUBROUTINE  FOR  RELATING  RTA  DATA  TO  SPL 
$FGRT 

C 

C THIS  SUBROUTINE  RELATES  THE  RTA  DATA  TO  A REFERENCE  SIGNAL 
C OF  KNOUN  AMPLITUDE. 


C 


C 

c 

c 

c 

c 

c 

c 

c 


FORTRAN  CALL  STATEMENT: 

-CALL  SPLCAL  (SPLCOR, NMAX,X, X0, NX0) 


SPLCOR  (REAL) 

NMAX  (I  NT) 

X (REAL) 

X0  ( INT*2) 

NX0  (I NT) 


CORRECTION  TO  BE  ADDED  TO  SUBSEQUENT  RTA  DATA 
TO  GET  SPL 

MAX.  NO.  OF  RTA  SAMPLES  ALLOWED  BY  CALLING  PROGRAM 
ARRAY  CONTAINING  RTA  FLOATING  POINT  DATA  FOR  1 BAND 
ARRAY  CONTAINING  THE  FIXED-POINT  RTA  DATA 
NO.  OF  ELEMENTS  IN  THE  ARRAY  X0; 

SHOULD  BE  CHOSEN  SO  THAT  NX0= (30*NSAMP) +60 


C 


10 


20 


30 

40 

50 

60 

70 

80 


90 


100 


OTHER  SUBROUTINES  CALLED:  SAMPLE, BAND. MEANLG 


SUBROUT I NE  SPLCAL ( SPLCOR . NMAX, X, X0 , MX0) 
DIMENSION  X(NMAX) 

INTEGERS  X0  (NX0) , TX100.  OVFLOLJ,  BANDNO 
REAL  MEAN 


URITE  (5,10) 

FORMAT  CTO  RELATE  THE  RTA  DATA  TO  A KNOUN  REFERENCE  SIGNAL,' 
l'TUO  CALIBRATION  OPTIONS  ARE  AVAILABLE : ' /X/ 

2'  OPTION  1:  USE  AN  ACOUSTIC  CALIBRATING  DEVICE  AS  A'/ 

3'  REFERENCE  SOUND  SOURCE' /X/ 


4'  OPTION  2:  MEASURE  THE  RMS  VOLTAGE  OF  A REFERENCE  SIGNAL, 
5'  COMPUTE  THE  SOUND  PRESSURE  LEVEL  FROM  THE'/ 

6'  SYSTEM  SENS ITI VITY'/X/ 

7 ' DO  YOU  UANT  TO  USE  OPTION  1?') 

READ  (5,20)  GS 


FORMAT  ( A4) 

IF  CGS.EQ. 'NO')  GO  TO  50 
URITE  (5,30) 

FORMAT  ('ENTER  THE  REFERENCE  SOUND  PRESSURE  LEVEL 
1 IN  DB  (XXX. XX)') 

READ  (5,40)  SPL 
FORMAT  CFG. 2) 

GO  TO  110 
URITE  (5,60) 

FORMAT  ('APPLY  CALIBRATION  SIGNAL  TO  SYSTEM,  AND  MEASURE  RMS 
1 VOLTAGE' /'OF  SIGNAL  AT  INPUT  TO  RTA  MULT  IF ILTER' /X/ 

2' ENTER  RMS  VALUE  IN  VOLTS  (X.XXXXX)') 

READ  (5,70)  VRMS 
FORMAT  ( F 7 . 5) 


URITE  (5,80) 


FORMAT  CENTER  SYSTEM  SENSITIVITY 
1 (-XX.XX)') 

READ  (5,90)  SENS 
F U R MA  T ( F 6 . 2 ) 


IN  DB  RE  1V/N/SQ. 


METER 


SPL =20*ALOG 10 (VRMS/. 00002) -SENS 
URITE  (5,  100) SPL 

FORMAT  (X/' SOUND  PRESSURE  LEVEL  CORRESPONDING  TO  CALIBRATION 
1 S I GNAL ' /' = ' , F7 . 2 , ' DB  RE  . 00002  N/SQ . METER ' /X/ 

2' DO  YOU  UANT  TO  CONTINUE?') 


READ  (5,20)  GS 
IF  CGS.EQ. 'NO')  GO  TO  5 


/ 


AND'/ 
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SPLCAL-2 


61  110 
62  120 


64 

65 

66 
67 
6S 

69 

70 


130 

135 

140 


7 1 


73  150 

74  160 

75 

7 b 

77  170 

78  175 

79 

80 
81 
DO 

'-JL- 

8 -• 

84  180 

85 

86 

87 

88 

89 

90  190 

91 

92  200 

93 

94 

95  210 

96 

97 

98 

99 
100 
101 

102  220 

103 

104 

105 

106  230 

107 

108 
109 
1 10 

111  240 
1 12 

113  250 
1 14 
1 15 
1 16 

117 

1 18  260 
1 19 

120  270 


URITE  (5,120) 

FORMAT  CENTER  THE  BAND  NO.  CORRESPONDING  TO  THE  FREQ. 

1 OF  THE  CALIBRATION  SIGNAL'/ 

2'  (IF  BAND  14  IS  USED  TO  SAMPLE  UNFILTERED  SIGNAL,  ENTER  "14")') 
READ  (5,130)  BANDNO 
FORMAT  (12) 

URITE  (5,140) 

FORMAT  ('HAVE  YOU  MADE  ANY  CHANGE  IN  THE  SYSTEM  GAIN  TO 
1 ACCOMMODATE' /'THE  CALIBRATION  SIGNAL?') 

READ  (5,20)  GS 
DGA  IN=0 . 

IF  (GS.EQ.'NO')  GO  TO  170 
LJR I TE  (5,1 60 ) 

FORMAT  CENTER  TOTAL  AMOUNT  OF  GAIN  CHANGE  IN  DB 
1 (MINUS  = REDUCTION)'/' ( -XX. XXX)') 

READ  (5,90)  DGA IN 
URITE  (5,175)  NMAX 

FORMAT  ('ENTER  TOTAL  RTA  SAMPLING  TIME  DESIRED  FOR 
1 REFERENCE  SIGNAL  MEASUREMENT'/ 

2' (IN  SECONDS,  2 DIGITS,  MAX  = ',12,')') 

READ  (5,130)  NS AMP 

IF  (NS AMP. GT. NMAX)  GO  TO  170 

URITE  (5,180) 

FORMAT  ('ENTER  "GO"  TO  START  SAMPLING') 

READ  (5,20)  GS 
TX100= 100 

CALL  SAMPLE (0, NSAMP,TX1 00, OVFLOU,X0,NX0) 

IF  (OVFLOU. EQ . 0)  GO  TO  200 
URITE  (5,190) 

FORMAT  ('READJUST  SYSTEM  GAIN') 

GO  TO  15U 

CALL  BAND (BANDNO, NSAMP, X, X0, NX0) 

CALL  MEANLG (NSAMP, X, MEAN, STDEV) 

URITE  (5,210)  MEAN, STDEV 

FORMAT  ('AVERAGED  RTA  READING  = ',F6.2,'  DB'/ 

1' STANDARD  DEVIATION  = ',F6.2,'  DB'/X/ 

2' DO  YOU  UANT  TO  REPEAT  THE  RTA  SAMPLING?') 

READ  (5,20)  GS 
IF  CGS.NE.'NO')  GO  TO  135 
FRECOR=0 . 

URITE  (5,220) 

FORMAT  ('DO  YOU  UANT  TO  CORRECT  FOR  FREQUENCY  RESPONSE?') 

READ  (5,20)  GS 
IF  (GS.EQ.'NO')  GO  TO  240 
URITE  (5,230) 

FORMAT  ('FROM  SYSTEM  FREQUENCY  RESPONSE  CALIBRATION  DATA, 

1 DETERMINE'/ 

2' RESPONSE  (CALIBRATION  FREQUENCY)  - RESPONSE  (1000  HZ)'/X/ 

3' ENTER  THIS  DIFFERENCE  (-XX.XX)') 

READ  (5,90) FRECOR 
SPLCOR=SPL-MEAN+FRECOR+DGA IN 
URITE  (5,250)  SPLCOR 

FORMAT  ('THE  CORRECTION  TO  BE  ADDED  TO  THE  RTA  READINGS  = ', 

IF 7. 2,'  DB'/X/' REMOVE  CALIBRATION  SIGNAL') 

IF  (DGAIN.EQ.0. ) GO  TO  270 
DGA IN=- 1 . *DGA IN 
URITE  (5,260)  DGA  IN 

FORMAT  C ADJUST  SYSTEM  GAIN  BY  ',F6.2,'  DB  (BACK  TO 
1 ORIGINAL  SETTING)') 

URITE  (5,280) 
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SPLCAL-3 


280  FORMAT  ('ENTER  "GO"  TO  CONTINUE*  OR  "STOP"  TO  REPEAT 
1 CALIBRATION') 

READ  (5*20)  GS 

IF  (GS.EQ. 'STOP' ) GO  TO  5 

RETURN 

END 
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FLOW  CHART:  SUBROUTINE  MEANLG 
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SUBROUTINE  MEAN LG 


$AS8M 

MEAN LG  PROG  MEANLG  - SUBROUTINE  FOR  FINDING  LOGARITHMIC  MEAN  AND  STD 
T-FORT 

C 

C THIS  SUBROUTINE  COMPUTES  THE  LOGARITHMIC  MEAN  AND 
C STANDARD  DEVIATION  OF  A NUMBER  OF  SAMPLES 
C FORTRAN  CALL  STATEMENT: 

C -CALL  MEANLG  (N, X, MEAN, STDEV) 

C N ( I NT)  = NUMBER  OF  SAMPLES 

C X (REAL)  = ARRAY  OF  SAMPLES 

C MEAN  (REAL)  = MEAN  VALUE  OF  THE  SAMPLES 

C STDEV  (REAL)  = STANDARD  DEVIATION  OF  THE  SAMPLES 

C 

C'  NO  OTHER  SUBROUTINES  CALLED 

C 

SUBROUT  I NE  MEANLG ( N , X, MEAN , STDEV) 

DIMENSION  X ( N ) 

PEAL  MEAN 
0=0. 

DO  10  1=1, N 
C=C+10.**(X( I)/10. ) 

10  CONTINUE 

U=C/N 

MEAN= 10. *ALOG 10(C) 

C =0 . 

DO  20  I = 1 , N 

C=C  + (X(  I)-MEAN)*CX(  D-MEAN) 

20  CONTINUE 

IF  (N.EG.l)  GO  TO  30 
C=C/(N- 1 ) 

30  STDEV=SQRT (C) 

RETURN 

END 
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FLOW  CHART:  SUBROUTINE  MEANAR 


C 


ENTER 


) 


V 

f 

COMPUTE 
1 N 
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\ 

r 
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COMPUTE 

_ 1 

' N 1 

WTJ  £ ^Xl 

1 ~ 1 -J 

\ 

f 

c 


RETURN 


) 
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SUBROUTINE  MEANAR 


1 $ASSM 

2 MEANAR  PROG  MEANAR  - SUBROUTINE  FOR  FINDING  ARITHMETIC  MEAN  AND  STD 

3 TFORT 

4 C 

5 C THIS  SUBROUTINE  COMPUTES  THE  ARITHMETIC  MEAN  AND 

6 C STANDARD  DEVIATION  OF  A NUMBER  OF  SAMPLES 

7 C FORTRAN  CALL  STATEMENT: 

8 C -CALL  MEANAR  (N,X, MEAN, STDEV) 

9 C N ( I NT)  = NUMBER  OF  SAMPLES 

10  C X (REAL)  = ARRAY  OF  SAMPLES 

11  C MEAN  (REAL)  = MEAN  VALUE  OF  THE  SAMPLES 

12  C STDEV  (REAL)  = STANDARD  DEVIATION  OF  THE  SAMPLES 

13  C 

14  C NO  OTHER  SUBROUTINES  CALLED 

15  C 

16  SUBROUTINE  MEANAR ( N , X, MEAN , STDEV) 

17  DIMENSION  XCN) 

18  REAL  MEAN 

19  C=0 . 

20  DO  10  I = 1 , N 

21  C-C+XCI) 

22  10  CONTINUE 

23  MEAN=C/N 

24  C = 0 . 

25  DO  20  1=1. H 

26  C=C+(X(  I)-MEAN)*CXU)-MEAN) 

27  20  CONTINUE 

28  IF  (N.EQ.l)  GO  TO  30 

29  C=C/(N-1) 

30  38  STDEV=SQRT (C) 

31  RETURN 

32  END 
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FEDERAL  INFORMATION  PROCESSING  STANDARD  SOFTWARE  SUMMARY 


01.  Summary  date 


Yr. 

Mo. 

Day 

1 1 _ 

1 

La, 

04.  Software  date 

Yr. 

Mo. 

Day 

u 


112. 


210. 


02.  Summary  prepared  by  (Name  and  Phone) 

Thomas  W.  Bartel  (301)  921-3783 


05.  Software  title 

Interactive  Program  for  the  Determination  of 
Sound  Power 


06.  Short  title 


SPOWER 


03.  Summary  action 

New  Replacement  Deletion 

ED  □ 

Previous  Internal  Software  ID 


07.  Internal  Software  ID 

None 


08.  Software  type 

09.  Processing  mode 

10. 

Application  are 

General 

Automated  Data 

Computer  Systems 

Management/ 

Q System 

[X]  Interactive 

Q Support/Utility 

□ Business 

[X]  Computer  Program 

Q Batch 

fc]  Scientific/Engineering 

□ Process  Control 

□ Subroutine/Module 

^ Combination 

[j  Blbllographlc/Textual 

□ Other 

Specific 

Acoustics;  reverberation 
room  measurements 


11.  Submitting  organization  and  address 


12.  Technical  contact(s)  and  phone 


Institute  for  Basic  Standards 
National  Bureau  of  Standards 
Washington,  D.C.  20234 


Thomas  W.  Bartel 
(301)  921-3783 


13.  Narrative 


The  program  is  used  to  measure  the  acoustic  sound  power  output  of  a source  in  a rever- 
beration room  according  to  ANSI  SI. 21-1972,  American  National  Standard  Methods  for  the 
Determination  of  Sound  Power  Levels  of  Small  Sources  in  Reverberation  Rooms.  The  com- 
puter controls  elements  of  the  instrumentation  system  through  a special  interface,  ac- 
quires the  digitized  sound  pressure  levels  of  the  reverberant  sound  field,  and  computes 
the  sound  power  levels.  Transmission  of  data  to  and  from  the  instrumentation  system  is 
handled  through  separate  subroutines  written  in  the  computer’s  assembly  language. 


14.  Keywords 

Acoustics;  computer-controlled  instrumentation  system;  real-time  analyzer; 
reverberation  room;  sound  power. 


15.  Computer  manuf'r  and  model 

16.  Computer  operating  system 

17.  Programing  language(s) 

Fortran  V,  Level  1 

18.  Number  of  source  program  state- 
ments 

Interdata  Model  70 

RT0S 

Interdata  Assembly  lang 

1058 

19.  Computer  memory  requirements 

20.  Tape  drives 

21.  Disk/Drum  units 

22.  Terminals 

37,066  8-bit  bytes 

None 

One/2.5  megabyte  unit 

One 

23.  Other  operational  requirements 


132-character  line  printer 

real-time  analyzer,  microphone  multiplexer  interfaced  to  computer 


1 t-.  V-  J-  ) uiiCiopuv 

24.  Software  availability 

25.  Documentation  availability 

Available 

s 

Limited 

□ 

In-house  only 
□ 

Available  Inadequate 

a □ 

In-house  only 
□ 

Paper  or 

mag  tape 

Internal  report 

26.  FOR  SUBMITTING  ORGANIZATION  USE 


Standard  Form  185 

1974  July 

U.S.  Dept,  of  Commerce — NBS 
(FIPS.  Pub.  30) 
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INSTRUCTIONS 


01.  Summary  Date.  Enter  date  summary  prepared.  Use  Year,  Month,  Day  format:  YYMMDD. 

02.  Summary  Prepared  By.  Enter  name  and  phone  number  (including  area  code)  of  individual  who  prepared  this  summary. 

03.  Summary  Action.  Mark  the  appropriate  box  for  new  summary,  replacement  summary  or  deletion  of  summary.  If  this  software  summary  is  a replace- 
ment, enter  under  “Previous  Internal  Software  ID"  the  internal  software  identification  as  reported  in  item  07  of  the  original  summary,  and  enter  the 
new  internal  software  identification  in  item  07  of  this  form;  complete  all  other  items  as  for  a new  summary.  If  a software  summary  is  to  be  deleted, 
enter  under  "Previous  Internal  Software  ID"  the  internal  software  identification  as  reported  in  item  07  of  the  original  summary;  complete  only  items 
01,  02,  03  and  11  on  this  form. 

04.  Software  Date.  Enter  date  software  was  completed  or  last  updated.  Use  Year,  Month,  Day  format:  YYMMDD. 

05.  Software  Title.  Make  title  as  descriptive  as  possible. 

06.  Short  Title.  (Optional)  Enter  commonly  used  abbreviation  or  acronym  which  identifies  the  software. 

07.  Internal  Software  ID.  Enter  a unique  identification  number  or  code. 

08.  Software  Type.  Mark  the  appropriate  box  for  an  Automated  Data  System  (set  of  computer  programs),  Computer  Program,  or  Subroutine/Module, 
whichever  best  describes  the  software. 

09.  Processing  Mode.  Mark  the  appropriate  box  for  an  interactive,  Batch,  or  Combination  mode,  whichever  best  describes  the  software. 

10.  Application  Area. 

General;  Mark  the  appropriate  box  which  best  describes  the  general  area  of  application  from  among: 

Computer  Systems  Support/Utility  Process  Control 

Management/Business  Bibliographic/Textual 

Scientific/Engineering  Other 

Specific:  Specify  the  sub-area  of  application;  e g.:  "COBOL  optimizer”  if  the  general  area  is  "Computer  Systems  Support/Utility";  “Payroll"  if  the 
general  area  is  "Management/Business";  etc.  Elaborate  here  if  the  general  area  is  "Other." 

11.  Submitting  Organization  and  Address.  Identify  the  organization  responsible  for  the  software  as  completely  as  possible,  to  the  Branch  or  Division  level, 
but  including  Agency,  Department  (Bureau/Administration),  Service,  Corporation,  Commission,  or  Council.  Fill  in  complete  mailing  address,  including 
mail  code,  street  address,  city,  state,  and  ZIP  code. 

12.  Technical  Contact(s)  and  Phone:  Enter  person(s)  or  office(s)  to  be  contacted  for  technical  information  on  subject  matter  and/or  operational  aspects 
of  software.  Include  telephone  area  code.  Provide  organization  name  and  mailing  address,  if  different  from  that  in  item  11. 

13.  Narrative.  Describe  concisely  the  problem  addressed  and  methods  of  solution.  Include  significant  factors  such  as  special  operating  system  modifi- 
cations, security  concerns,  relationships  to  other  software,  input  and  output  media,  virtual  memory  requirements,  and  unique  hardware  features. 
Cite  references,  if  appropriate. 

14.  Keywords.  List  significant  words  or  phrases  which  reflect  the  functions,  applications  and  features  of  the  software.  Separate  entries  with  semicolons. 

15.  Computer  Manufacturer  and  Model.  Identify  mainframe  computer(s)  on  which  software  is  operational. 

16.  Computer  Operating  System.  Enter  name,  number,  and  release  under  which  software  is  operating.  Identify  enhancements  in  the  Narrative  (item  13). 

17.  Programing  Language(s).  Identify  the  language(s)  in  which  the  software  is  written,  including  version;  e.g.,  ANSI  COBOL,  FORTRAN  V,  SIMSCRIPT  11.5, 
SLEUTH  II. 

18.  Number  of  Source  Program  Statements.  Include  statements  in  this  software,  separate  macros,  called  subroutines,  etc. 

19.  Computer  Memory  Requirements.  Enter  minimum  internal  memory  necessary  to  execute  software,  exclusive  of  memory  required  for  the  operating 
system.  Specify  words,  bytes,  characters,  etc.,  and  number  of  bits  p>er  unit.  Identify  virtual  memory  requirements  in  the  Narrative  (item  13). 

20.  Tape  Drives.  Identify  number  needed  to  operate  software.  Specify,  if  critical,  manufacturer,  model,  tracks,  recording  density,  etc. 

21.  Disk/Drum  Units.  Identify  number  and  size  (in  same  units  as  "Memory" — item  19)  needed  to  operate  software.  Specify,  if  critical,  manufacturer, 
model,  etc. 

22.  Terminals.  Identify  number  of  terminals  required.  Specify,  if  critical,  type,  speed,  character  set,  screen/line  size,  etc. 

23  Other  Operational  Requirements.  Identify  peripheral  devices,  support  software,  or  related  equipment  not  indicated  above,  e.g.,  optical  character  devices, 
facsimile,  computer-output  microfilm,  graphic  plotters. 

24.  Software  Availability.  Mark  the  appropriate  box  which  best  describes  the  software  availability  from  among:  Available  to  the  Public,  Limited  Availability 
(e.g.:  for  government  use  only),  and  For-ln-house  Use  Only.  If  the  software  is  "Available",  include  a mail  or  phone  contact  [joint,  as  well  as  the  price 
and  form  in  which  the  software  is  available,  if  possible. 

25.  Documentation  Availability.  Mark  the  appropriate  box  which  best  describes  the  documentation  availability  from  among:  Available  to  the  Public,  In- 
adequate for  Distribution,  and  For  In-house  Use  Only.  If  documentation  is  "Available",  include  a mail  or  phone  contact  point,  as  well  as  the  price  and 
form  in  which  the  documentation  is  available,  if  possible.  If  documentation  is  presently  "Inadequate”,  show  the  expected  availability  date. 

26.  For  Submitting  Organization  Use.  This  area  is  provided  for  the  use  of  the  organization  submitting  this  summary.  It  may  contain  any  information 
deemed  useful  for  internal  operation. 
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